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(54) Abstract Title 

Angle dependent surface multiple attenuation for two • 



component marine bottom sensor data 



(57) An angle dependent filter for two-component 
sensor data allows for attenuation of free surface 
multiples. Typical sensors are hydrophones and 
geophones. The method decomposes the recorded dual 
sensor data into upgoing and downgoing wavefields by 
combining the recorded pressure at the hydrophone with 
the vertical particle velocity from the geophone recorded 
at the ocean floor. Surface multiple attenuation is 
accomplished by application of an incident angle 
dependent deconvolution of the downgoing wavefield 
from the upgoing wavefield. The method uses an angle 
dependent filter to calibrate the geophone response so 
that the different coupling of the two instruments and 
associated noise are taken into account. 

In a further embodiment, a method of attenuation 
of multiple reflections in seismic data is provided. The 
seismic data comprises pressure data and particle 
velocity data. The method comprises deconvolving the 
seismic data, and applying a moving average operator to 
the seismic data. 
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Transform the 
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(for each p trace) 
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each p trace) 
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ANGLE DEPENDENT SURFACE MULTIPLE ATTENUATION FOR 
TWO-COMPONENT MARINE BOTTOM SENSOR DATA 



The present invention relates generally to marine seismic exploration, and more 
particularly to a marine seismic measurement system that allows for attenuation of free 

10 surface multiples in two-component marine bottom sensor data. 

Marine seismic wave measurement systems are used to take seismic profiles of 
underwater geological configurations. One procedure of marine seismic measurement 
involves the use of a marine bottom cable. Surveys using marine bottom cables are often 
employed in areas that are populated with numerous obstacles, such as drilling and 

15 production platforms. In a procedure using marine bottom cable, bottom cables are 
deployed along the marine bottom. Often, multiple cables are deployed in parallel. Each 
bottom cable has a plurality of sensor pairs placed at regular intervals along the cable. 
Each sensor pair contains a pressure sensor, such as a hydrophone, and a particle velocity 
sensor, such as a geophonc. A gimbal mechanism within each geophone ensures that the 

20 sensing elements of the geophones arc vertically oriented. 

Acoustic energy is generated in the vicinity of the marine bottom cables using an 
acoustic energy source such as an air gun array or a marine vibrator array. An air gun 
discharges air under very high pressure into the water. Marine vibrators typically 
include a pneumatic or hydraulic actuator that causes an acoustic piston to vibrate at a 

25 range of selected frequencies. The vibrations of the acoustic vibrator produce pressure 
differentials in the water that generate acoustical energy pulses. Source acoustical waves 
travel downward through the water and into the earth as seismic waves. The, source 
waves strike interfaces between formations in the earth. A portion of the source wave is , 
reflected upwards from the interface towards the marine bottom. The sensor array on the 

30 marine bottom receives the reflected waves and converts the waves into signals that are 
recorded as sensor data. The sensor data is processed to provide information about the 
structure of the formations beneath the marine bottom. 

The sensor array receives not only the reflected waves of interest, but also the 
source waves and reverberated waves. Reverberated waves are waves that have been 



reflected from the water-air interface back towards the marine bottom. Such 
reverberated waves are referred to as free surface multiples or surface multiples. The 
free surface multiples may be significant in amplitude and may be difficult to 
differentiate from the desired reflected waves. 

The use of dual sensor measurements, namely pressure and vertical particle 
velocity, allow for the attenuation of free surface multiples. U.S. Pat. Nos. 5,163,028; 
5 365 492- 5 524,100, and 5,621,700 describe methods of attenuating free surface 
m ultiples,'and each of these patents are incorporated by reference as if fully sot forth 
herein The methods of attenuating free surface multiples detailed in the above 
0 referenced patents do not adequately take into consideration the angle dependence of the 
upgoing and downgoing wavefields. Also, the methods of attenuating multiples detailed 
in the above referenced patents do not adequately take into consideration the angle 
dependency of the response of each sensor of a sensor pair. The use of methods of 
attenuating multiples that do not consider both the angle dependency of the upgoing and 
,5 downgoing wavefields and the angle dependency of the response of each sensor of a 
sensor pair may lead to inaccurate determinations of the formations present beneath the 
marine bottom. 

Existing methods of deconvolution for multiple attenuation of dual sensor data 
carry out the calibration and deconvolution filter in the distance-time domain, (x.f). The 
20 basic equations for deconvolution of upgoing and downgoing waves are valid only m the 
angle (plane wave) domain. Seismograms recorded from a single shot will have energy 
propagating at all possible angles, so processing data in the distance-time domain can 
only have limited success. 

SUMMARY OF THE INVENTION 
The problems outlined above are in large part to be solved by a system and 
method of marine seismic exploration that takes into account angle dependencies dunng 
the processing of two-component sensor data. Consideration of the angle dependences 
of upgoing wavefields, downgoing wavefields, and the sensor enhance the attenuate of 
free surface multiples that are present in two-component sensor data. The ability to 
30 provide enhanced attenuation of free surface multiples allow for more accurate 
determination of the formations present beneath a marine bottom. 

In one example embodiment, a method is used to decompose the recorded dual 
sensor data into upgoing and downing wavefields in the plane wave domain. The 



mothod finds an angle dependent calibration factor that allows the calibration of the 
recorded pressure data with respect to the recorded vertical particle velocity data. The 
angle dependent calibration factor takes into consideration the angle dependencies of the 
hydrophone and the geophone, as well as noise associated with the recording geometry. 
Attenuation of multiples is accomplished by application of an incident angle dependent 
deconvohrtion of the downgoing wavefield from the upgoing wavefield calculated from 
the calibrated pressure and vertical particle velocity data. 

In a further embodiment, a method of attenuation of multiple reflections in 
seismic data is provided. The aeismic data comprises pressure data and particle velocity 

data. The method comprises deconvolving the seismic data, and applying a moving 

average operator to the seismic data. 

In a further embodiment, a system for attenuation of multiple reflections in 

seismic data is provided. The seismic data comprises pressure data and particle velocity 

data. The system comprises means for deconvolving the seismic data, and means for 

applying a moving average operator to the seismic data. 

In a further embodiment, seismic data is provided. The seismic data is processed 

by a method comprising deconvolving the seismic data, and applying a moving average 

operator to the seismic data. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Further advantages of the present invention will become apparent to those skilled 
in the art with the benefit of the following detailed description of example embodiments 
and upon reference to the accompanying drawings in which: 

FIG. 1 is an illustration of a marine seismic survey system; 

FIG. 2 is an illustration of a primary wave impinging a sensor pair; 

FIG. 3 is a flow diagram representing a method of processing marine seismic 
survey data that takes into consideration the angle dependency of the sensor response and 
the angle dependency of the upgoing and downgoing waveftelds; 

FIG. 4a and 4b show pressure and vertical velocity data, reapectively, from 
synthetic shot data for a 1-D four layer model; 

FIG. 5a and 5b show the data represented in Figures 4a and 4b transformed to the 

plane wave domain; 

FIG. 6a and 6b show the calculated separated upgoing wavefield and the 
corresponding calculated downgoing wavefield, respectively. 



FIG. 7a shows the multiple attenuation results in the distance-time domain; 

FIG 7b shows the location of the simulated primaries of the model; 

FIG. 8 shows hydrophone data from an offshore ocean bottom experiment in the 

South China Sea; 

FIG. 9 shows corresponding geophone data for the ocean bottom experiment in 
the South China Sea depicted in Figure 8; and 

FIG. 1 0 shows the multiple attenuation results obtained after processing the data 

represented in Figures 8 and 9. 

DETAILED DESCRIPTION OF EXAMPLE EMBODIMENTS OF THE PRESENT 

INVENTION 

With reference to the drawings, end particularly to Figure 1, a uuu-ine aeiamic 
surV ey system is generally designated by reference numeral 10. In one embodiment of 
the present invention, the system 10 includce a seismic survey ship 12 una, tows an 
acoustic energy source 14 through a body of water .6. In various embedments, the 
rustic energy source is an array of acoustic energy sources. In alternate embodiments, 
the acoustic energy source is an air gun, a marine vibrator, or another dev.ee that 
generates acoustic waves. The construction and operation of acoustical energy .sources ■ 
wel! knmvn in the art and is no. described in detail herein. The activation of the acousuc 
enercy source is refelTed to as "shooting." 

,n a further embodiment, the system 10 includes receiving ship 18. The receivtng 
ship 18 deploys bottom cable 20 on marine bottom 22. The receiving ship deploys an 
amy of bottom cables in parallel lines. Each bottom cable 20 carries a. least one sensor 
pair 24, and preferably, each bottom cable carries a plurality of sensor pairs. Each sensor 
pair 24 includes a pressure sensing transducer, such as a hydrophone, and a parade 
velocity sensor, such as a geophone. As is well known in the art, each marine geophone 
indudes a gimbal mechanism to ensure tat the sensing dement of the geophone ts 
vertical., oriented during use. In a fitrther embodiment, each hydrophone and geophone. 
sends separate data signab to Ore receiving ship IS. The data is recorded by. mult, 
channe, seismic recording system tha, selective* amplifies, condthons. and record, 
. rime-varying electrical signals. In sill, a further embodiment, the system also digmzes 
„« received signals to facilitate signal analysis. Any of a variety of aeiamic recordmg 
systems is used to record the data. 
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To take a marine seismic survey, the receiver ship 18 positions the bottom cable 
20 on the marine bottom 22. In an embodiment, shooting takes place while the survey 
ship 12 moves at a constant speed along a set of survey lines with respect to the cable 20. 
The location and depth of each sensor pair 24, and the location of the acoustical energy 
source 14 at the time of each shot are recorded. After the survey ship 12 completes the 
survey line, the receiving ship 1 8 retrieves the cable 20 and redeploys the cable in a new 
position. After re-deployment of the cable 20, the survey ship 12 shoots another set of 

survey lines. 

During data collection, seismic waves 26 generated by the source 14 travel away 
from the source. Portions of the waves travel downward and into the land beneath the 
marine bottom 22. The waves are reflected off of interfaces between subterranean 
formations, such as interface 28 between subterranean formations 30 and 32 as shown in 
Figure 1. Reflected waves 34 from the interfaces travel upwards and impinge upon a 
sensor pair 24. The sensor pairs 24 detect the reflected waves 34 and transmit signals 
along the cable 20 to the receiving ship 18. The receiving ship 18 records the data so 
that the data can be subsequently processed to map the location of interfaces 28 between 

subterranean formations. 

The sensor pairs 24 receive not only the reflected waves 34, which are also 
known as primaries, but also the source waves 26 and free surface multiples 36. The free 
surface multiples 36 may be significant in amplitude and may be difficult to differentiate 
from the desired reflected waves 34. Free surface multiples that originate from the 
source, contact the air/surface interface 38, and travel towards the marine bottom are 

referred to as ghosts 40. 

In most marine seismic data acquisition situations, the energy source 14 is placed 
above the marine bottom 22. All upgoing wavefields in the data result from reflections 
of downgoing incident waves to the marine bottom. Mathematically, this is formulated 
by the following convolution, 

d(tr = r(t)®d(t) dn (D 

where d(t) up are the upgoing wavefields, which can be either primaries 34 or multiples 
36; and d(t) d " is the downgoing wavefield, which may be direct transmission from the 
source 26 or reflection events 36 (source, primaries or multiples) that are bounced back 
at the air/surface interface 38 (receiver ghosts). In equation (1), r(t) is the reflectivity of 



.he shuCure. which includes .hose reflections taking p.ace both at the manne bottom 22 
1 inside ,h= .end — beiow ,hc marine bottom: and . is «* wavcfidd h—e. 
The reflectivity, r(t). also includes internal multiple retlecfons. 

The reflectivity, r, is obtained by deconvo.ving the downgomg wavefield from 
; the upgoing wavefield. In one example embodiment, in ,he frequency domain, tins is 
represented by a simple division: 

R(a.)-D(<u^/D(a>)'" , (2) 
where R(»), WO", and DO* ere the Fourier transforms with respect ,o time of r«). 

HfT, and d(t)"\ respectively. 

The deconvolution based multiple attenuation method given m eouatton (2) 
quires separated upgoing and doling wavefiOds. The upgoing and downgomg 
wLie.ds are repre.en.ed by the upgoing and downgoing pressure, or by the upgomg 
and downgoing vertical particle velocity. 

Z the above description, the angle dependency of the propagating waves was 
15 suppressed to present the fundamen..! ideas behind processing two-componen. manne 
bottom sensor data. The deve.npmen, below Wees into account angie dependence. 

In various embodiments, seismic data is recorded in tire offset d.etance-.une 
domain. Offset distance is the hor.aon.ai distance between the lection of a sensor pa,r 
^ end .he .oeerien of the eeoostic energy source M at U.O time Of a shot. The S«Snu 
20 data recorded b, a sensor pair 24 is arecord of tite variation in pressure, as measured by 
Mrophones, and vertica. partic.e ve.ocity, as measured by geophones mWen as 
function of sourcc-to-receiver offset distance and time. In a tether embodnnent, * 
^stic energy source ,4 is considered to be a point source, m practice, me so.ee M 
wiB have a directivity that is angie dependent In stil. a further 
2S response of a source is synthesized by summing the responses from a ser.es of ptane 
JL each eharac.er.ed by the propagation angle of the piane ^ In 
embodiment, the source generated dara recorded in tire offeet-trme doma-nr 

data where * is the vector mpresenting the source-.o-receiver offset distance, and . - me 
30 frequency, the plane wave response is given by: 



u(T,p3 = Jdffi.dxu(a>,e) 



where E is the vector ray-parameter, and t is the offset time, or x = t-p. x. In a two- 
dimensional geometry, p=p x =p=sin9/a. The angle 9 is the angle of propagation, and a is 
the velocity of sound in the medium. The angle of propagation 9 is illustrated in Figure 2 
for a primary wave 34. 

S In a further embodiment, in the frequency-ray-parameter domain, the pressure 

recorded by hydrophones is given by the equation: 

P(co,p) - C 1 (P)(1+R(^.P))- 1 + R(".P>- iu,,,h )S(co,p). (3) 

The corresponding vertical velocity geophone data in the frequency-ray parameter 
I o domain is given by the equation: 

V 2 «o,p) - c 2 ( P )(KR( ffi ,p))' , (R(",Pr (,h - e- to « h )S(a>,p), (4) 

where p is the ray-parameter. 9 is the angle of propagation; a is the sound velocity in the 
medium; R is the reflectivity of the structure referenced to the water surface, which 

15 includes the response of both the primaries and internal multiples; h r is the receiver 
depth; q is the vertical slowness, or q = «l/a W); S is the source wavelet with source 
side ghosts; and c, and c 2 are two incident angle dependent coefficients defined by the 
ray-parameter and the parameters of the medium. 

In a further embodiment, the upgoing and downgoing wavefields are obtained by 

20 combining the pressure and the vertical components of the vertical particle velocities in 
equations (3) and (4) to yield: 

P(p)"P « (1/2) (P(p) + c(p)V a (p)), P(p) dn = (1/2) (P(P) - c(p) V 2 (p)) , (5) 

where c(p) is a function defined by the parameters of the medium. 

25 In some embodiments, equation (5) is considered to be an equation describing a 

system wherein the hydrophones and the gcophones are perfectly coupled to the 
environment, wherein the geophones have the same response characteristics, and wherein 
the environment surrounding the sensor pair 24 is noise free. Direct application of 
equation (5) may not yield good results because, in reality, a sensor pair is not perfectly 

30 coupled to the environment, the instrument response characteristics of hydrophones and 
geophones are not the same, and sensor pairs 24 are not located in noise free 
environments. In further embodiments, the angle dependency of the sensor pairs require 
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that the geophone data be calibrated to the pressure data, or vice versa, before the 
upgoing and downgoing wave components are determined. 

In a further embodiment, from equations (3) and (4), the following equation are 

derived: 

S(co,p) = (l/2)[(l+Z(cD,p))3P(a>,p)[(l-Z(a>,p))]V 2 (a,,p) (6) 

where Z(co,p) is a filter and §(co,p) is the temporal Fourier transform of the time delayed 
source function, which is a function of the source excitation function S(a>,p), and is grven 
by the equation: 

S(oa,p) = (l/C2qa J ))e'^ h (a** - e** S(co,p)) 
The term §(a>,p) contains the source excitation function S«o.p), the source side ghosts, 

and the transmission operator. 

Since the aource excitation function is always of finite duration in real data, there 

must exist a time T 0 (p), such that: 

s(t,p) = 0,.ift>t 0 + h 1 /a < 7 ) 

where s(t,p) is the inverse temporal Fourier transform of §(*#). which is the delayed 
source ftinction* 

In a further embodiment, equation (7) allows the hydrophone data and the 
geophone data to be calibrated so that the delayed source function is optimal. The 
20 calibration filter is designed such that the delayed source function defined in equation (7) 
will have minimum energy after a certain time. The time will include the source 
excitation function time duration and the sum of the propagation time of the energy of 
the source to go to the receiver and back to the surface. In a further embodiment, the 
calibration function,.^), is found by solving the following equation and constraint: ^ 

IJF'Kl+Z^^PCto^+Xto.pXl-Z^'P))!^^!! =minimum (8) 
t>T 0 = to+h r /a 

where F» stands for the inverse Fourier transform operator. In a further embodiment, a 
numerical method of solving for the calibration function in a least squares sense involves 
30 solving a system of equations with Toepliz structure to find a value for the calibrate 
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function for a given value of p. A person having ordinary skill in the art will recognize 
that in alternate embodiments, many methods are used to solve the equation and 
constraint of equation (8) to yield the optimum angle dependent calibration function. 

Finding the calibration function allows for the calibration of hydrophone and 
geophone data. In yet a further embodiment, angle dependent multiple attenuated data is 
then obtained by forming and using an optimal filter Prom the calibrated data by 
application of equations (5) and (2). 

The methodology for one embodiment of processing two-component marine 
sensor data is shown diagrammatically in Figur e 3. Figure 3 shows that the data is first 
collected and recorded in the offset distance-time domain. The data is then transformed 
to the plane wave domain, x- P . Each p trace of the transformed data set is converted into 
multiple free data by the application of the steps illustrated within block 42. For each p 
trace of the transformed data set, a calibration function is numerically calculated based 
on equation (8). The calibration function is then used to calibrate the vertical particle 
velocity data with respect to the pressure data. The calibrated data is then used to find 
information corresponding to the upgoing and downgoing wavefields. The upgoing and 
downgoing wavefield information is used to numerically solve for an optimal filter to 
suppress multiples by applying Equation (2). The optimal filter for the P trace is then 
used to generate multiple free data for the p trace. The generated multiple free data for 
each P trace forms a multiple free data set. The multiple free data sets are transformed 

back to the offset time domain. 

Figures 4-7 show an application of one example embodiment of processing two- 
component marine seismic data to a 1-D four layer acoustic model. The hydrophone and 
geophone data are shown in Figures 4a and 4b. Figure 5a and 5b show the t-p 
transformation of the data. The pressure shown in Figure 5b is then decomposed into an 
angle dependent estimate of the upgoing wavefield and the downgoing wavefield by 
applying equation (5). The upgoing wavefield is shown in Figure 6a, and the downgoing 
wavefield is shown in Figure 6b. The reflectivity of the model is computed by 
deconvolving the downgoing wavefield from the upgoing wavefield in the x-p domain. 
Following an inverse t-p transform, the final result in the x-t domain is obtained and 
shown in Figure 7a. By comparing the results shown in Figure 7a with the original data 
sets shown in Figures 4a and 4b, and with the simulated primary reflections of the model 
as shown in Figure 7b, it is seen that the method of processing two-component marine 
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^ data attenuates tire multiples in the data and recovers the reflect of tt. 

,„ *«« «. *. a,™ «, *. - *. *~ — » 

which repr^en, the watcr-iand interface, and two subterranean interfaces. 

One exampie embodimen, of the method has aiso been apphed to ocean bottom 
, doa! sensor dam coUected in me South China Sea. Figure 8 shows me P~>^^ 
Figore 9 shows the corresponding vertical particte vetoed After transforming ,o *. 
p,L wave domain, the gephone data is eaUbrated before aeperanag the up^ and 
downgoing wavef.eUs. The final rest* of nudtiple attenuation rs shown m Figure 10. 
I IpLon wrtb the origin,, dam, the resuh has been convolved with . 
. wave,,, derived from me autocorrelation of each p .race in the hydrophone da, before 
Lrse e-p hansfom, As shown in Figure 10, me muhip.es have been — * - 
Lore importantly, the primary reflections, which are marked b, arrows, have been 

PreSCtV t primary reflections were no, singly visibie in ,he origmal data due ,o me 

15 presence of strong multiples. ...... 

,„ an even further embodimen, of me presentinvention, me s.gna, - eond., oaed 
hefore me multiple eUmination procedure is applied. In an alternate embodtmen,, me 
sigaaHs conditioned after the multiple elimination procedure is apphed 

„ stiti a further embodimen, ang.e dependence is contented ,» tirevertrcal 
20 «M* P-e wave seismogram and me pressure piane wave seismogram. 

,. ffo „ nr(M , ; n waV e types recorded and m the recoraing 
in the response are due to the differences in wave ty P 

SMK ors. As described above, in various embodiment minimizing *e source^ 
dotation defines an angle dependen, filter, which is used te cahbrate tire fwo p.ane wave 

2S a further embodiment a convention, predictive deconvomtion is apphed 

,o ,he plane wave pressure and vertica, seismograms. ,n a firrther 
predictive deconvomtion is applied in me plane wave domain. In alternate embody 
L deconvCution is applied te tire date .ogemer or independentiy. «n one : emb^menU 
the read, is ma, each will now have tire short time duration component of 
30 dependent responses, which arise from each sensor (geophone and hydrophone) 

"""t still a further embodimen,, foUowing me angle dependen, deconvom.io.cach 
seismogram is bandpass fihered ,o tire same frequency band. In ye, a further 
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embodiment, an amplitude calibration is performed for the deconvolved and bandpass 
filtered plane wave pressure and vertical velocity seismograms. In one embodiment, the 
amplitude calibration uses time signal moving average operators for each plane wave 
component to define the mean signal level. In further example embodiments, the time 
5 signal moving averaging is achieved through the following equations: 

Ik 7='-* 

In still a further embodiment, an inter-seismogram calibration further comprises 
10 comparing an average of pressure seismograms with an average of vertical velocity 
seismograms for each pair (pressure and vertical velocity) of plane wave seismograms. 
In yet a further embodiment, equations (9) and (10) are used in equation (8) for 
calibration. 

In another embodiment, the vertical velocity signal average is normalized to that 
15 of the pressure average to define the amplitude correction on a time sample by time 
sample basis. In a further embodiment, the time signal moving average operator length k 
is varied to achieve the optimum calibration. Predictive deconvolution followed by 
moving averaging of signals applied to plane wave seismograms introduces stability and 
results in robust estimates of multiple free seismograms. 
20 In a further embodiment, a method of attenuation of multiple reflections in 

seismic data is provided. The seismic data comprises pressure data and particle velocity 
data. The method comprises deconvolving the seismic data, and applying a moving 
average operator to the seismic data. 

In a further embodiment, said deconvolving the seismic data further comprises 
25 applying a predictive deconvolution to the seismic data. In a further embodiment, h 
applying a predictive deconvolution to the seismic data further comprises applying a 
predictive deconvolution to the pressure data and the particle velocity data 
independently. In alternate embodiments, a deconvolution is applied such as a spiking 
deconvolution, gapped deconvolution. statistical deconvolution. deterministic 



12 



dcconvolulion or any other devolution .hat will occur to those of ordinary s W U in the 

^ In still a further embodiment, the method further comprises bandpass fihenngthe 
seismic data. In a further embodiment, .he method mrther comprises bandpass filtenng 
, the oressure data and <hc particle velocity data independently. 

Ta tamer ombodimen, me method comprises amplitude calibrating me ae^ue 
da*, In a further embodiment, said amplitude ea.ibro.ing mrther comprises defimng a 

""tlH. embcdimenr, said applying a moving average operator further 
,0 comprises applying essentially the following formula: 

ta a further embodiment, said applying a moviae average operator mrther 
, 5 comprises applying essentially me following formula: 

„ still a further embodiment, a system for at.enua.ion of multiple reflections in 

• • a ♦ ^TT^rkes nressure data and particle velocity 
seismic data is provided. The seismic data composes pressure data p 

a,„« for deconvolving the seismic data, and means for 
™ data The system compnses means for deconvolving 

2£, moving average operator to the seismic da,, in — embodiment. Ore 
m dec.nvo.ving .he seismic data comprises a compute, a worsen softwar. 
. sof ,ware running on any computer, or any o.her means m, dcconvo.vmg .be «-» da* 
th a, wm occur to those of ordinary Ml in the art. In aliemate embodiments. *™ 

orator to the seismic data comprises a computer, a 
?<i for aDolvmc a moving average operator to me sc. 

; I .L software, software running on any compute, or any Cher mean s * 
applving a moving average operamr to .be seismic data .ha, w... occur to those 

ordinary skill in the art. r**u~r 
b a further embodiment, said means for deconvolving me sersnuc da* fnrth r 
,o comprises means for applying a predictive deconvohrtion to the seismic data In 
Za.« enrbodimen., the means for app.,ing a predictive deconvo.ution .o ft. se.sm.c 
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data comprises a computer, a workstation, software, software running on any computer, 
or any other means for applying a predictive deconvolution to the seismic data that will 
occur to those of ordinary skill in the art. 

In a further embodiment, said means for applying a predictive deconvolution to 

5 the seismic data further comprises means for applying a predictive deconvolution to the 
pressure data and the particle velocity data independently. In alternate embodiments, the 
means for applying a predictive deconvolution to the pressure data and the particle 
velocity data independently comprises a computer, a workstation, software, software 
running on any computer, or any other means for applying a predictive deconvolution to 

10 the pressure data and the particle velocity data independently that will occur to those of 

ordinary skill in the art 

In a further embodiment, the system further comprises means for bandpass 
filtering the seismic data. In alternate embodiments, the means for bandpass filtering the 
seismic data comprises a computer, a workstation, software, software running on any 
15 computer, or any other means for bandpass filtering the seismic data that will occur to 

those of ordinary skill in the art. 

In a further embodiment, the system further comprises means for amplitude 
calibrating the seismic data. In alternate embodiments, the means for amplitude 
calibrating comprise a computer, a workstation, software, software running on any 
20 computer, or any other means for amplitude calibrating that will occur to those of 

ordinary skill in the art. 

In a further embodiment, said means for applying a moving average operator 
further comprises means for applying a time signal moving average operator. In 
alternate embodiments, the means for applying a time signal moving average operator 
25 comprises a computer, a workstation, software, software running on any computer, or 
any other means for applying a time signal moving average operator that will occur to 

those of ordinary skill in the art. 

In a further embodiment, said means for applying a moving average operator to 
the seismic data further comprises means for inter-sei sinogram averaging. In alternate 
30 embodiments, the means for inter-seisraogram averaging comprises a computer, a 
workstation, software, software running on any computer, or any other means for inter- 
seismogram averaging that will occur to those of ordinary skill in the art. 
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In a further embodiment, means for inter-seismogram calibration further 
comprises means for comparing an average of the particle velocity data and the pressure 
data In alternate embodiments, the means for comparing an average of the particle 
velocity data and the pressure data comprises a computer, a workstation, software, 
5 software running on any computer, or any other means for comparing an average ol the 
particle velocity data and the prcasure data that will occur to those of ordinary sktll » the 
art. 

In a further embodiment, the system further comprises means for normalizing the 
particle velocity data. In alternate embodiments, the means for normalizing comprises a 
,o computer, a workstation, software, software running on any computer, or any other 
means for normalizing that will occur to those of ordinary skill in the art. 

In a further embodiment, the system further comprises means for normalizing the 
pressure data. In alternate embodiments, the means for normalizing comprises a 
computer, a workstation, software, software running on any computer, or any other 
, 5 means for normalizing that will occur to those of ordinary skill in the an. 

in a further embodiment, seismic data is provided. The seismic data is processed 
by a method comprising deconvolving the seismic data; and applying a moving average 
operator to the seismic data. 
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CLAIMS 

1 . A method for attenuation of multiples in marine seismic data, comprising: 

creating seismic wavefields in a water environment, each created wavefield 
5 generated at a different location; 

recording sensor readings from at least one sensor pair, wherein said at least one 

sensor pair comprises two different types of sensors, wherein said sensor readings are 

substantially a result of the created seismic wavefields, and wherein said sensor reading 

form a data set; 

10 calculating angle dependent estimated values for upgoing and downgoing 

wavefields that generated the data set; and 

generating an attenuated multiple data set using the estimated values for the 
upgoing and downgoing wavefields. 

15 2. The method as defined in claim 1, wherein at least one sensor pair comprises a 
pressure sensor and a vertical particle velocity sensor. 

3 The n etted of claim 1 or claim 2, vterein calculating angle dependent estinotBd 
values for upgoing and downgoing wavefields further comprises: 
20 transforming the data set to a plane wave domain data set; 

calibrating the plane wave domain data set using angle dependent calibration 
functions; and 

using the calibrated plane wave domain set to calculate the estimated values for 
the angle dependent upgoing and downgoing wavefields. 

25 

4. A system for generating attenuated multiple data for two-component marine 
seismic sensor data, comprising: 

an acoustic energy source configured to generate an acoustic wavefield; 

at least one sensor pair, the at least one sensor pair comprising two types of 
30 sensors configured to detect wavefields resulting from the generation of an acoustic 
wavefield by the acoustic energy source; 

a recording system configured to separately record sensor readings from each 
sensor of the at least one sensor pair as data, the recording system configured to record 
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th e lo ca«ion of ,he acoustic energy source a, a time u*en *e acoustic : energy source 
generates art acoustic wavefeld, and said recording system configured to record the 
location of the at least one sensor pair, and 

. data processor configured to anaiyze the data and produce the attenuated 

multiple data; , ~ 

wherein the data processor calculates angle dependent estimated values of 

upgoing and downgoing wavers that generated the data set, and wherein the 
processor uses the estimated values of the upgoing and downgoing wavefields to 
generate the attenuated multiple data during use. 

5. A method of generating attenuated multiple data from two-component marine 

seismic data, comprising: 

i j estimated values for upgoing wavefields and 

calculating angle dependent estimated vam * y& 

downgoing wavefields that generated the data; and 

generating an attenuated multiple data set using the estimated values-for the 

upgoing and downgoing wavefields. 

6 A method of attenuation of multiple reflections in seismic data, wherein the 
seismic data comprises pressure data and particle velocity data, the method compnsing: 

, deconvolve the seismic data in the plane wave domain; and 

applying a moving average operator to the seismic data. 

7 The method of claim 6, wherein said deconvolving the seismic data further 
comprises applying a predictive devolution to the seismic data. 

, The method of claim 7, wherein said applying a predictive devolution to the 
seismic data further comprises applying a predictive devolution to the pressure data 
and the particle velocity data independently. 
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9 . TrenHthodof any of clirins 6 to 8. further oarprising bandpass filtering tie 
seisnic data. 

,0. The method of claim 9. further comprising bandpass filtering the pressure data 
and the particle velocity data independently. 
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] ] t The method of any c£ claims 6 to 10, aarprising arplitude <3librating the 
seisrric data, 

12. The method of claim 11, wherein said amplitude calibrating further comprises 
defining a mean signal level. 

5 

13 The method of any of claims 6 to 12, \sherein said applying a nrving average eperatrn 
further adiposes applying a time signal noting average c^eidluL using essentially the 
fhl lowing fbmulae: 

2K J=i-k 

14. The netted of aiy of claims 6 to 13, vterein said ^plying a roving average qperatar 
to the seismic data further emprises imer-seiarogran calibration 

15 15. The method of claim 14, wherein said inter-seismogram calibration further 
comprises comparing an average of the particle velocity data and the pressure data. 

16. The nethod af any of claims 6 to 15, further crnpcising nanializing the 
particle \elccity and pressure cfeta* 

20 

17. A system for attenuation of multiple reflections in seismic data, wherein the 
seismic data comprises pressure data and particle velocity data, the system comprising: 

means for deconvolving the seismic data in the plane wave domain; and 
means for applying a moving average operator to the seismic data. 

25 

18. The system of claim 17, where said means for deconvolving the seismic data 
further comprises means for applying a predictive deconvolution to the seismic data. 
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1*. The systen of claim 17 car claim 18, further oaTprising maans for besx^ass 
filtering the seismic data. 
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20. He system of any of claims 17 to 19. corprising means fear aipLLtude 
calibrating tte ssisrac date. 

21 . The system of claim 20, wherein said means for amplitude calibrating furdier 
5 comprises means for defining a mean signal level. 

22 He system of cry of claims 17 to 21, therein said neans fear applying a ncving average 
operator to the seismic data further comprises means for inter-seismogram calibrate 

l0 23 The system of claim 22, wherein said means for inter-seismogram calibration 
. * mT ,«rtnP an average of the particle velocity data and the 

further comprises means for comparing an average u* f 

pressure data. 

24. The system of any of claims 17 to 23, further cotprising maans tor roaLLzing toa 
j 5 particle velocity and pressure data. 

25 Seismic data processed by a method comprising: 

deconvolving the seismic data in the plane wave domain; and 
applying a moving average operator to the seismic data. 

20 
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